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Renormalisation of Newton’s constant 


Kevin Falls 
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Philosophenweg 12, 69120 Heidelberg, Germany 

The problem of obtaining a gauge independent beta function for Newton’s constant is addressed. 
By a specific parameterisation of metric fluctuations a gauge independent functional integral is 
constructed for the semiclassical theory around an arbitrary Einstein space. The effective action 
then has the property that only physical polarisations of the graviton contribute, while all other 
modes cancel with the functional measure. We are then able to compute a gauge independent 
beta function for Newton’s constant in d-dimensions to one-loop order. No Landau pole is present 
provided N g < 18, where N g = d(d - 3)/2 is the number of polarisations of the graviton. While 
adding a large number of matter fields can change this picture, the absence of a pole persists for 
the particle content of the standard model in four spacetime dimensions. 
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I. INTRODUCTION 

In quantum field theory it is well known that coupling constants become functions of the energy scales 
entering the renormalisation process. In turn this implies a modification of the classical scaling properties 
of a theory [l|. Such energy dependence of a coupling a is encoded in its beta functions f3 a = where /r 

is the renormalisation scale. In quantum chromodynamics (QCD) this non-trivial scaling can be observed 
in the one-loop beta function [2j, J3j], 

«-"( u -!*')£■ (11) 

where a s is the strong coupling and Nf denotes the number of flavours. For Nf <16 this equation describes 
the weakening of the strong force as the energy scale n is increased. On the other hand if we have Nf > 16 the 
coupling constant will diverge at a finite energy. In the former case the theory is said to be asymptotically 
free and is well defined for all scales. In the later case the coupling has a Landau pole and the theory can 
only be considered as an effective one, at energies below the pole. 

When considering gravity as a quantum field theory it is natural to ask whether there exists an analog 
to m for Newton’s constant Gn. In four spacetime dimensions Newton’s constant is dimensionful which 
implies that the structure of the beta function differs from that of QCD, where a s is dimensionless. In 
particular, for d spacetime dimensions one should work with the dimensionless coupling G = n d ~ 2 GN- The 
analog of m will then take the form 


/3 c = (d-2)G-6G 2 , (1.2) 

for some constant b. Whereas the coefficient of G just reflects the dimensionality of Gat the coefficient b 
is a quantum correction h Clearly b plays the same role as 11 - |iV/ in the QCD; if b is positive then 
G will remain finite for all scales, but if b is negative G will blow up at a finite energy. Thus, although 
non-perturbative effects may change the naive conclusion, the sign of b is a strong indiction of whether 
quantum gravity should be replaced by some new theory at a finite energy scale, or whether gravity may 
itself constitute a perfectly valid quantum field theory. The latter scenario, dubbed asymptotic safety by 
Weinberg [dj , would involve a fixed point for G where /3g = 0 for some positive value G = G*. At the level of 
the one-loop beta function (11.211 this is found at G* = j;(d- 2) and its existence only depends on the sign of b. 
However, one must go beyond one-loop in order to be sure that such a fixed point exists non-perturbatively. 

Although /T; has been calculated previously, in these calculations the value of b has been found to depend 
on the gauge [5j-l8( and parameterisation [8j, |9( (an exception are calculations with the geometrical effective 
action [lOj] , which are gauge independent although a particular parameterisation of the fields is used). Hence 
these calculations will generically include unphysical degrees of freedom coming from the gauge fixing action 
and may not be trust worthy. This issue stems from the fact that calculations have to be performed off-shell 
and therefore generally depend on both the gauge and the parameterisation. When working on-shell such 
dependencies vanish. On the other hand the Einstein equations dictate that the scalar curvature is given 
by 

(1.3, 

where A is the cosmological constant. As a consequence one cannot disentangle the renormalisation of 
Newton’s constant from the vacuum energy, forcing us to work off-shell to determine b. Thus, it is clear 
that by working off-shell, in a gauge dependent setup, no definitive conclusion can be reached as to the sign 
of b. On the other hand explicit non-perturbative calculations both in four and higher dimensions indicate 


1 b is proportional to h. although here we work in natural units h = 1 = c. 
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that a fixed point for G exists SB [M3, despite calculations being generally gauge dependent. Hence, 
although evidence is strong that the coupling G does not blow up for some finite energy scale, it is also 
questionable since the dependence on the gauge parameters implies that unphysical contributions are still 
present. 

In this paper we shall obtain gauge independent results by disentangling physical degrees of freedom 
at the level of the functional integral. In particular we shall obtain a semiclassical approximation to the 
d-dimensional functional integral in quantum gravity on an arbitrary Einstein space independent of the 
gauge without fixing R by the on-shell condition m- This functional integral has the property that it 
only receives contributions from physical fluctuations of the metric that survive on-shell. This leads us to 
gauge independent beta function for the gravitational constant. We will then see that the sign of b can be 
universally determined by a unique factor 


b°c 1(18-N g ), 


with 


N g = 


d(d - 3) 


(1.4) 


where N g is the number of dynamical degrees of the metric in general relativity. Thus one finds that b is 
positive in d = 4 dimensions but becomes negative for integer dimension d > 8. This result is for pure gravity 
but is also easily generalised when matter is included. 

At this point we should take a moment to clarify the meaning of (11.41) . Although one may like to think 
of G being related to an observable, this may not be the case. In particular, at least perturbatively, there 
is no universal definition of a running Newton’s constant for e.g. scattering processes |l8[. What ED does 
tell us is whether or not the ultra-violet (UV) cutoff A can be removed and hence whether a continuum 
limit exists. For N g < 18 one can remove A -»■ oo without any problems. However, for N g > 18 one may show 
that for the bare Newton’s constant to be positive, and hence for the functional integral to exist, A must 
remain finite. Thus the beta function implies the existence (non-existence) of a continuum limit for N g < 18 
(N g > 18). The universal and gauge independent result (11.41) therefore provides evidence for the continuum 
limit of quantum gravity in four spacetime dimensions. 

The rest of this paper is as follows. In section El we find the form of the functional measure for the 
gauge fixed functional integral over geometries and use it to determine the on-shell functional integral for 
the semi-classical theory. This generalises the result of 0 to all dimensions d > 2. Turning to section ED 
we then confront the problem of obtaining gauge independent results off-shell finding a parameterisation 
of metric fluctuations which achieve this end given in section IIII Al We are then able to write down a 
gauge independent one-loop effective action. In section IIII Bl the gauge independence is linked to the gauge 
invariance of the quadratic action achieved while using our specific parameterisation. Since our gauge 
independent result is achieved using a transverse-traceless decomposition in section IIII Cl we also calculate 
the effective action in de Donder gauge as a cross-check and find a form of the effective action in terms 
of unconstrained fields. Section m is devoted to a calculation of quantum corrections to the Newtonian 
potential showing that the result is independent of the parameterisation. We then use our parameterisation 
to derive a semiclassical renormalisation group equation for a scale dependent action T u in section El In 
section El we evaluate traces appearing in the flow equation using the early time heat kernel expansion and 
access their universal content. We may then observe the universal factor (11.41) while also reproducing the 
curvature squared counter term found in [2^ • Then in section IVIII we give the explicit form of the beta 
function /3g■ We first discuss the case near two dimensions and the limit d -*■ 2. Then we consider the 
general case in d dimensions and discuss in more detail the physical origin of the universal factor (11.41) . In 
section rVIIIl we find the form of the divergent counter term for the vacuum energy and note that it vanishes 
in d = 3 dimensions. In section m we derive a bound on the UV scale A of the theory in higher dimensions 
and show no bound exists for N g < 18. In Section |X] we present the beta function with the inclusion of 
matter. We end with our conclusions in section m 
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II. THE ON-SHELL FUNCTIONAL INTEGRAL 


General relativity in d dimensions describes (d- 3)d/2 dynamical degrees of freedom coming from the 
d(d+ l)/2 components of the metric, which are determined by Einstein equations, minus d constraints and 
d diffeomorphisms. In the quantum theory these degrees of freedom correspond to the N g polarisations of 
the graviton, schematically one has 


N g = d(d+ l)/2 


d = d(d - 3)/2. 


metric diffeomorphisms constraints 


( 2 . 1 ) 


Thus one should expect that the additional 2d unphysical degrees of freedom are removed when physical 
quantities are computed. At the level of the functional integral such cancelations occur between metric 
fluctuations and the functional measure. One way to ensure this cancelation occurs is to turn off any exter¬ 
nal source terms which violate parameterisation invariance, equivalently, when working with the effective 
action, this can be achieved by working on-shell. In this section we will determine the form of the semi- 
classical functional integral by expanding around a solution to the equations of motion to ensure gauge 
and parameterisation independence. This we do as an intermediate step to gain physical intuition and see 
mathematically how the unphysical states are removed when computing on-shell quantities. 

The first step in our calculation is to determine the form of the functional measure for quantum gravity. 
We assume a bare action for Eunclidean quantum gravity of the Einstein-Hilbert form 

5srav[0>] = 16^ / d d XsJ 7[2A b - R (7 

flV ) ] + Sgfbv], (2-2) 


where S g { is the gauge fixing action. Here 7 ;ii , is the metric tensor with determinant 7 = det 7 /i „ and R{ 7 ^) 
denotes the Ricci scalar. The action depends on two parameters Gb and A b which denote the bare Newton’s 
constant and bare cosmological constant respectively. The functional integral is then given by 

Z = J I? 7 AlJ/ (detQ)e“ Sgrav[7 ' 1 " ] (2.3) 


where det Q is the determinant of the Faddeev-Popov operator. 

Here we shall evaluate the functional integral (12.311 in the semiclassical approximation. To this end we 
first consider parameterisations of metric fluctuations h^ in terms of a background field g^ which satisfies 
the classical equations of motion 2 


= d-2 


A b 9[iv ’ 


Two such parameterisations are the linear parameterisation 


(2.4) 


lliv - 


[J/ju h/j 


(2.5) 


and the exponential parameterisation 


"/fj.u - <7/iA[e ]j/ 



( 2 . 6 ) 


where in the later case h is a symmetric matrix with components \h\£ 5 h£. An advantage of the exponential 
parameterisation (12.61) is that one may single out the conformal factor of the metric by 


= ei 9^l eh ]t = 


(2.7) 


2 Here we adopt the notation that all curvatures and covariant derivatives are with respect to unless indicated otherwise, 
such that R = R(g^ u ) etc. and V/x = V/x(p7zu) is the covariant derivative with respect to g^u, and indices are lowered and 
raised with g^u- 
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where <f> - h'p is the trace of h and h is the trace free part. From here it follows that the determinant of 7 
is fixed to the background one 


7 = 9 (2-8) 

<f> 

and thus is identified as the conformal factor of 7 ^„ with respect to the background metric g^ v . 

We shall need to expand the action to second order in the fluctuation around the on-shell metric g^. 
To ensure the resulting hessians are invertible we must also fix a gauge, the choice of which should fall out 
of physical quantities. Here we choose a class of background field gauges, linear in the fluctuation for 
either parameterisation, 

% = 32^ / ddx V9d^ F ^, ^ = Va^-^V^. (2.9) 

We now follow the steps out lined in [2l| to decompose the fluctuation such that it becomes manifest 
that the gauge fixing action only depends on d independent fields matching the number of diffeomorphisms. 
The first step is to adopt the transverse-traceless decomposition of the fluctuations in terms of differentially 
constrained fields , 

h/iv — + 4 > ~^g^i' ^ 1 ( 2 . 10 ) 

h^ = 0 , = 0 . = 0 . 

One observes that takes the form of a transverse diffeomorphism of the metric to linear order. Thus we 
can identify ^ as d - 1 of the unphysical fields corresponding to such diffeomorphisms. We then further 
redefine the scalar fields {cp,^} -* {</>, ip} as 

for which it becomes manifest that the gauge fixing action (12.91) only depends on the transverse vector 
and the scalar ip. Thus we see that ip represents the additional longitudinal diffeomorphism. One should 
then expect that the integral of these fields in (12.31) should be cancelled by the functional measure. The 
field redefinition (12.101) leads to the following Jacobians in the functional measure 

Jh — y/ deto[A]y/deto [Ao]\/det it [Ai], (2-12) 

whereas the redefinition (12.111) has a trivial Jacobian. The non-trivial Jacobian (12.121) is given by determi¬ 
nants of the differential operators, 


A ip=-V 2 ip 

A i'p li = (-V%-R lt v )<p v , 

acting on scalars and vectors as indicated. The subscripts IT and 0 specify that the determinants are 
evaluated from transverse vectors and scalars. It is important to note that the constant mode of the scalar 
ip should be left out of the functional measure since it cannot contribute to the fluctuation h Accordingly 
the the constant modes in the scalar Jacobians should also be removed. Similarly one should leave out the 
Killing vectors = 0 from the vectors and the transverse vector Jacobian ^/det [Ai], this 

corresponds to the zero mode of Ai. 

Next we turn to the Faddeev-Popov determinant. Here we will exploit the liberty to write the determinant 
as det Q = Vdet Q 2 such that the ghosts become second order in derivatives at the price of including a set 
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of commutative ghosts fields B^ in addition to the anticommuting ghosts and We then decompose 
the ghosts into transverse and longitudinal fields 

= Bj + = Cj + , C p = Cj + V M C L , (2.14) 

where the superscript T indicates that the field is transverse i.e. V^Bj = 0. This transformation then leads 
to a Jacobian J g h = /1 which cancels with the corresponding factor in the metric fluctuation Jacobian 

Jh given by (12.1211 leaving just two Jacobians, 

Jo = \/deto [Ao], J\ = \J detir [Ai]. (2-15) 

It is these Jacobians that remove the d degrees of freedom corresponding to the constraints whereas the 
ghosts should cancel the integral over diffeomorphisms (i.e. and ip). Due to the decomposition (12.141) 
the corresponding Faddeev-Popov determinants splits into two factors for the longitudinal and transverse 
parts, 

det Q L - yj det 0 [A|A], det Q T = ^/detrr[^?] (2-16) 

where we define the gauge dependent differential operator by 

(2 - 17) 

acting on scalar fields. The path integral is then given by 

Z= J Vh^VpV^Vip J 0 J^etiQ^detiQr) . (2.18) 

To obtain the semiclassical functional integral we must then expand the action to second order in the 
fluctuations of the fields {<p, while taking into account the contributions from the measure. The 

next stage is therefore to compute the second variation of the bare action and apply the on-shell condition 

(E9). 

First let’s show that the choice of parameterisation, such as m or (USD, does not effect the on-shell 
hessians. Here we just assume that the metric is a function of the fluctuation h A (x) (with e.g. A - {y^}) 
such that 


7 /J.V = hny{h A {x)) with 7 m „(0) = g^ 


(2.19) 


which agree only up to linear order in the fluctuation. The choices (12.61) and (12.51) are just two examples of 
such parameterisations. Introducing de Witt’s condensed notation we write i = {x, gv} and a = {x, A } and 
take two functional derivatives with respect to the fluctuations h A (x) = h a expressed with derivatives with 
respect to the metric "/^(x) = 7 *. Acting on e.g. <S'[ 7 I ] = *S gr av we obtain, 

j 2 B[ 7 fc [fe a ]] s fsy wn . ss[f] fry* sstf] 

Sh a Sh b 5h b \5h a S'y 1 ) Sh a h b Sj 1 Sh a Sh b S'y’S'yi 


on-shell the first term is zero since it is proportional to the equation of motion, while the second term only 
depends on 7 ^ to linear order of in h This shows that by going on-shell the results cannot depend on 
which of parameterisations is used. 

/ Q \ 

The on-shell hessians 5g r i v has the following components, for metric fluctuation fields ip and 

167rG & 4 2) = iA 2 , 167rG b 5g= (d ~J d ; p)2 A|A, 16 ttG 6 = ^A 2 , (2.21) 


5 2 S l 


grav 


5<p A (x)ip B (y) 


where 


VHS ( Xb lABS(x-y) 


( 2 . 22 ) 



















with ip A denoting the various fields and 1 ab is the identity on the corresponding field space. In (12.211) the 
differential operator A 2 is the Lichnerowicz Laplacian acting on two tensors as 

A 2 ^ = -VV-2H/ v P <p a p. (2.23) 

The field <j>, whose on-shell hessian is given by 

16 *G b S™ = , (2.24) 

needs extra attention since one observes that (12.241) is negative for positive eigenvalues of the operator A 0 . 
This indicates that the naive Wick rotation of the functional integral is unbounded from below. However, 
as pointed out by Mottola and Mazur [l9| one should Wick rotate all modes with positive eigenvalue as 
<p -*■ icj) whereas modes with negative eigenvalues should be Wick rotated trivially. This rule can be derived 
by considering the super-metric on the space of fluctuations h M! ,. Hence one can write the corresponding 
operator as 

16 *G b S™ = ( ^~ 2rf^~ 2> |A°I ’ (2 - 25) 

to get a well defined Euclidean functional integral. 

Expanding the action to second order in the fluctuation, and comparing S^ and given in (12.211) to 
(12.161) . one observers that the functional integrals over £ and ip cancel with the Faddeev-Popov determinants 
det Ql and det Qt■ All gauge dependence has therefore cancelled out. Furthermore comparing (12.251) with 
Jo we see that all modes of <j> apart from the constant mode 


= 0, 


(2.26) 


are cancelled by the Jacobian Jo given by (12.15|) (here we assume that the constant mode is the only mode 
for which A 0 has a negative eigenvalue, which is true at least for the d-sphere). This follows since, as we 
argued, the constant mode must be left out of Jo- The only contributions that remain are the (d-2)(l + d)/2 
transverse-traceless fluctuations, the Jacobian J 1; which comprises d -1 negative degrees of freedom, and the 
constant mode <pQ. Thus one is left with (d- 3)d/2 local degrees of freedom corresponding to the graviton 
and one global degree of freedom corresponding to a constant rescaling of the metric. The functional integral 
then reduces to the form 


Z = 


J dtpo J T>/i 1 v / deti T [Ai]exp -SgravDs] - .^\ Gb f d<lx V9 ^ </>oA 0 ^o j 


(2.27) 

this result is valid in all dimensions d > 2 and for all parameterisations (12.191) of the metric fluctuations. 
The remaining determinant can be expressed in terms of integral over auxiliary fields 

1 


Ji = V d etiT [Ai] = J T>C^T>c IJ ,T>c ll exp ~ ^ G J d d x^/g (c m Aic m + C^Ai£ m ) 


yw 




(2.28) 


where are anticommuting and is commuting and all are transverse. Here 5a 


action that combines with the gravity action in (12.271) to form a modified bare action S = S gra .v t ^ a ux 


is the auxiliary 
+ Saux- The 

semiclassical functional integral (12.271) generalises the d = 4 result of [19]. We note that, while the cancellation 
of the d diffeomorphisms and one of the constrained degrees of freedom is explicit, the cancellation between 
the remaining d - 1 unphysical transverse traceless fluctuations is implicit and implies that non-trivial 
cancelations between these modes and the Jacobian Ji should occur. It is clear therefore that one should 
not leave out the Jacobians from the functional, as was done in [23j], doing so implicitly includes d unphysical 
degrees of freedom and fails to reproduce the one-loop functional integral (12.271) . The importance of including 
such Jacobians to obtain the correct covariant measure has been stressed in [24] for the case of QED. 
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III. THE OFF-SHELL EFFECTIVE ACTION AND GAUGE INDEPENDENCE 

It is well known that the off-shell effective action in quantum gravity generally depends on both the 
gauge and the field parameterisation (see e.g. [25j] ). On the other hand since we observe which modes 
survive in the on-shell semiclassical approximation, and which modes cancel in the functional measure, we 
understand which fields are carrying physics, and which fields are only present due to us necessarily breaking 
diffeomorphism invariance. 

In order not to be forced to expand around a solution to the bare equations of motion we will now drop 
the condition (12.41) and instead simply require that the background is an Einstein space, 

R )J .v = ^g t i V R- (3.1) 

As pointed out in the introduction this is crucial if we are to be able to extract the renormalisation of 
Newton’s constant. 


A. Gauge independence 


Our aim is to preserve the gauge independence of the on-shell functional integral (12.271) such that physically 
meaningful results may be derived off-shell. The problem that we face is that by working off-shell the 
cancellations between the measure and the 2d unphysical degrees of freedom is not guaranteed. For a 
generic parameterisation terms proportional to the equation of motion arise in the hessians and in the 
presence of such terms the cancellations observed in the previous section will not occur. In particular the 
hessians have the general form, 

S^ 2 ) = S (2) + x(r~-^ A h ), (3.2) 


for any arbitrary parameterisation, where X is a matrix in field space that depends on the parameterisation 
given by 


d-2dS^ 
~2d~ d\ b 


(3.3) 


which then implicitly defines S^ 2 \ While X is dependent on the parameterisation S^ is parameterisation 
independent. 

At this point we make two observations, if X + 0 the gauge fixing dependence will remain and furthermore 
will possess negative eigenmodes depending on A b- Thus, only by picking a parametrisation for which 
X = 0 will the gauge fixing dependence vanish and the Gaussian integrals converge (assuming positivity of 
S for the chosen background). By picking a parameterisation for which XtO the corresponding effective 
action will not formally exist for certain values of Ab where negative eigenvalues of S^ are present, as well 
as being gauge dependent. We are therefore led to the conclusion that a parameterisation for which X = 0 
is the key to obtaining a physical meaningful gauge independent one-loop result. In turn this removes any 
dependence of the hessians on the cosmological constant. 3 

To find such a parametrisation we first note that by using the exponential parameterisation (12.61) the 
equation of motion will only appear in the hessians for if> and (f>. In fact, if we revert to the the parameteri¬ 
sation in terms of ?/> and <^>, (see (12.111) ') the equation of motion will only appear in the hessian for <f> due to 


3 We note that the dependence on the cosmological constant can also be removed by use of specific gauge fixing conditions |8] 
whereby the conformal fluctuations cj) are constrained. 
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(12.81) . this field being identified as the conformal factor of the metric. Then we may choose to parameterise 
this scalar degree of freedom in terms of a field er related to the conformal factor by 


± 

e d 




2 -d 
4 d 2 


a + ... 


(3.4) 


such that 



— 9h/us + 


9/US 



(3.5) 


where h is a traceless fluctuation which may be decomposed as (12.101) with </> = 0. It then follows that 
X = 0 since the volume element is linear in the field and hence no term proportional to the cosmological 
constant can arise in the Hessian. 

It should be noted that the transformation between <f> and a is singular at the point a = —2 and at this 
point vanishes. Thus the full non-perturbative treatment of the functional measure would require extra 
care involving the determinant of the field space metric to ensure that the measure is re-parameterisation 
invariant. In addition one should decide over which domain the functional integral should be performed 
and take care of possible singularities arising for a given parameterisation, such as at a = -2. Here we are 
concerned only with the semiclassical theory based on the expansion in small a to second order for which 
we can neglect the Jacobian which results from (13.41) . One should also stress the choice (13.41) is not unique. 
All that we require for X to vanish is that v /y has no term quadratic in the field 


X = 0 


S 2 

Sh^(x)Sh pr7 (y) 



(3.6) 


For any such choice the same semi-classical functional integral will be obtained. Additionally it may be 
natural to have all higher derivatives of the spacetime volume V = f d d x^/ 7 to vanish as is the case for (13.51) . 
In this case we have that all terms proportional to the equation of motion are absent in all higher vertex 
functions S^ for n > 2 . 

Adopting (13.51) we then perform the transformations ( 12 . 111 ) but now for er instead of (j>, 


a = <t> + V 2 ip, 


= il> + 


(d- 1 - p)V 2 + R 


(3.7) 


After this transformation the hessians then reduce to the form ( 12 . 211 ) and (12.251) without use of the equations 
of motion. One may expand around an arbitrary Einstein background 63 and again obtain a functional 
integral given by (12.2711 , the upshot being that the scalar curvature R is no longer fixed by the cosmological 
constant. 

The one-loop effective action T can then be obtained from (12.271) in the standard way by adding a source 
term and performing a Legendre transform. We then set ( h pv ) = 0 to obtain a functional of a single metric 
field <? M „. On an arbitrary Einstein space the gauge independent action T then takes the simple form 


r [g^\-S[g pv \ = ^STr[log <S (2) ] = ^Tr^log A 2 ] - ^TrirflogAi] + ^log^-^J?, (3.8) 

where S is the total bare action appearing in the functional integral, with all determinants expressed in terms 
of auxiliary fields and ghosts, and (STr) Tr denotes the (super)-trace. The subscripts 2 T 2 and IT denote 
that the traces are over transverse-traceless tensor modes and transverse vector modes, respectively, and 
the final term is the contribution of the constant mode of 4>. The simplicity of this effective action compared 
other gauge fixed actions is quite appealing, in particular it is independent of any scalar modes (apart from 
the constant mode </> 0 ). Most importantly (13.81) has the property that all unphysical polarisations of the 
graviton, that cancel on-shell, also cancel off-shell. 
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B. Gauge invariance 


One may guess that the gauge independence which is observed when using the parameterisation (13.41) is 
linked to an underlying gauge invariance. In turn the lack of gauge independence for other parameterisations 
such as (12.51) or m is due to the absence of gauge invariance. In fact this can be observed by looking at 
the quadratic action 

S qu , d = h-S {2 \g^)-h (3.9) 

for 423D or & where we use a condensed notation. Apply a gauge transformation 

“I" + Vi. (3.10) 


One can then show that S'quad is only gauge invariant if the metric v is chosen to be on-shell. This 
observation has been made by Deser and Henneaux [26j] for the linear parameterisation . However if one 
instead takes the parameterisation (13.51) we will be lead to a quadratic action of the form 


“S'quad — 2 ^'^ ^ ) ' V 5 


(3.11) 


with ip = {a, h^}. 


One may then show that this action is invariant under the linear gauge transformations 


+ ,Vq£ 


a —» o’ + 2 V a £° 


(3.12) 


for an arbitrary Einstein space m■ Thus the gauge independence of the one loop effective action (13.81) 
can be traced precisely to the gauge invariance of the underlying quadratic action (13.111) resulting from the 
parametrisation (13.51) . 


C. de Donder gauge and the effective action for unconstrained fields 


We have shown that a gauge independent one-loop effective action can be found at least for the gauges 
& However, the understanding that this gauge independence is due to an underlying gauge invariance 
gives confidence that the effective action should be the same for other gauge choices outside the two param¬ 
eter family (12.91) . To show this gauge independence we have utilised the transverse-traceless decomposition 
(12.101) which also leads to the Jacobians (12.151) . For consistency we now derive the effective action without 
use of the decomposition (12.101) or the corresponding decomposition for the ghosts (12.141) while working 
explicitly in the de Donder gauge for which we have p = | -1 and a = 1 in (12.91) . This gauge choice has been 
studied many times for the linear parameterisation and also recently for the exponential parameterisation 
( 9 ]. Here we will use the parameterisation (13.51) . The hessians for the metric fluctuations then take the form 


siV = 


d -2 
~c T 


-V - -R 


d -2 
4 d 


|A S 


(3.13) 


S!f - ii 2 (3.14) 

where we have Wick rotated cr -*■ ia to give the positive sign for all modes apart from the constant 
mode. The corresponding ghost operator is given by 

det Q = detiAi. 


(3.15) 
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Here we note that the neither the metric fluctuations nor the ghosts have any differential constraints. One 
then obtains the following one-loop effective action 

r[flv] - S[g^ v ] = |ti -T 2 log A 2 + ^Tr 0 log |A S | - Tri log Ai (3.16) 

where the first trace Tr 2 t is over traceless tensor fluctuations, the second trace Tro is over scalar fluctuations 
and the last trace Tri is over vector modes. This effective action has a more familiar form of metric 
fluctuations minus ghosts fluctuations. Although its form is different from that of (13.81) one can use the 
relation between traces of differentially constrained fields and of non-differentially constrained fields to show 
that (13.161) and (13.81) are equivalent for d > 2. Explicitly one has that the trace over transverse traceless 
modes is given by a trace over traceless modes minus a trace over vector modes 

Tr 2 T 2 /(A 2 ) = TY 2 r /(A 2 ) -TW(Ar) (3.17) 


and that a trace over transverse vector modes is given by a trace over vectors minus a trace over scalars 

ThT/(Ai) = Tr 1 f(A 1 ) - Tr'/(A s ). (3.18) 

Here the scalar trace has the constant mode removed. These relations (gj, 27] hold on a generic Einstein 
space where we assume that there are no (conformal) killing vectors. One can then use (13.171) and (13.181) 
to show that (13.161) is indeed equal to (13.81) . To bring the effective action into it’s simplest form we can use 
the relation between traceless tensor fluctuations and an unconstrained symmetric tensor 


Tr 2 /(A 2 ) = Tr 2 T /(A 2 ) + Tr 0 /(|A S |) 


(3.19) 


where Tr 2 is the trace over symmetric two tensor modes with eigenvalue of constant mode sign flipped. 
Finally we get the expression 

r [g^u] - Slg^] = ^Tr 2 log A 2 - Tiq log Ai (3.20) 

which has the intuitive form of a trace over symmetric tensor fluctuations minus the ghost fluctuations with 
a relative factor of two in front of the ghosts without any further constraints for either field. An alternative 
route to (13.201) can be found by introducing a fluctuation h at the beginning as 

h/ns — h/iv ~ ; & — h p // 5 (3.21) 

such that we have just a single symmetric fluctuation field h^. Then one finds the hessian in de Donder 
gauge is given by 

S$ = \c- A 2 ( 3 . 22 ) 

where the contribution to the effective action from the local operator C pvXp = \{g PlX g p,p + g* lp g >lv - g pv g pX ) 
can be absorbed into the functional measure. 

The effective actions (13.201) . (13.161) and (13.8p are all equivalent on a generic Einstein space. Care must be 
taken to treat Killing vectors and conformal Killing vectors if a background metric is chosen where they are 
present. In this paper we will side set this subtlety by assuming that no such modes are present. 


IV. INTERLUDE: PARAMETERISATION INDEPENDENCE OF THE NEWTONIAN 

POTENTIAL 

The main aim of this paper is to compute beta functions which are independent of the gauge fixing em¬ 
ployed. This we take as a prerequisite for the running couplings to be considered as physically meaningful. 
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To find the beta function for Newton’s constant we will look at the local divergencies of the effective action, 
however it is the non-local terms in the effective action which can be related to low scale observables in 
quantum gravity mm. When computing an observable the dependence on both the gauge and parame- 
terisation should fall out of the result even though a dependence will generally be observed at intermediate 
stages of the calculation. 

To check that the observables are indeed parameterisation independent we will calculate the corrections 
to the Newtonian potential felt by a (classical) test particle via corrections to its geodesic equation [3(| . 
generalising the result to an arbitrary parameterisation. These corrections come from two sources, namely 
the quantum correction to the Einstein equations and the quantum corrections to the geodesic equation itself. 
Only when both of these are taken into account does the dependence on the gauge and parameterisation 
cancel out. Other approaches to defining a physical Newtonian potential have been considered such as [31 1 
based on scattering potential, however the explicit gauge independence has not been shown. 

Here we consider a classical source of mass M and a test particle of mass m coupled to gravity with a 
bare action of the form 


R 


S = J d 4 x^y - M J yj-^^dy^dy'' -m J Z~1^dz> l dz v 


167T Gb 


(4.1) 


where it is assumed that to « M. The classical field Einstein equations therefore involve a non-zero energy 
momentum tensor T+ T000. After computing the one loop effective action quantum corrections to 
the Newtonian potential V are found via the equation of motion for z in the non-relativistic limit 


VV = 


d 2 z 

Id? 


(4.2) 


where the particle of mass M is taken to be a stationary source. 
In 


30.] the potential V was calculated for a one parameter family of gauges for which p = 1 and a 


is the free parameter. There the linear parameterisation was used whereas here we wish to analysis the 
parameterisation dependence for a fixed gauge. To this end we consider the four parameter family of metrics 
given 


Tm 1 ' — 0^ ^Xu + 02 h X X + 03 g^h Xp h\ p + 04 g^h x xh p , 


(4.3) 


while fixing the gauge to p = 1 and a - 1. The parameterisation m corresponds to the choice 0\ = 1, 
02 = ~ 2 , 03 = 0 , 04 = Here we repeat the work of [30| recalling only the main results and stating the 
differences arising from taking 0t + 0. 

To find the corrections to the Newtonian potential we need to calculated the non-local terms in the effective 
action. First the quantum corrections due to the heavy particle and gravity can be found neglectin g th e 
term proportional to to. With 0i = 0 these terms take the form r nl (/3j = 0) = rgl(/3j = 0) + T^{0i = 0) 13(1 


r^C/3* = o) = ■ 


967T 2 


J d^Xx/^dj 


^^log(-n)i?^ + ^R\og(-a)R 


(4.4) 


n\(0i = 0) = -g-^ f d 4 Xx/^[Mp Vpa \og(-n)M p ^ + 2Mp Vpa \og(-n)^P p ^ + ^Rg^g^ . 

where the tensors P and M are given by 

PZ = 2 Rx 25<X) - Z'Rx, - gxaR^ - RS^xK) + \z v 9x«R (4.5) 

MZ = M 2 4ttG n J drS 4 (x-y(T))[g IJ -''yxya + ‘2y IJ 'y''yxya] (4.6) 

Here we have calculated the non-local terms for general parameterisation finding r nl = r nl (/3j = 0) + 

= J d A XxZZWr)R pv \og{-a)£ pv + B(0 i )Rg^log(-a)£^] (4.7) 
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with 

A Pi G n e G n { 439ft + 58& + 1768& + 232fo) . . 

32tt ’ 384tt ' 1 ' ’ 

As expected the effective action for different parameterisations only differs by terms proportional to the 
field equations £ M1/ = - 167rGjy {R' lv - \Rg + \T^ (here we neglect terms quadratic in the equation of 
motion). In fact the term T^ 1 has the same form as terms which arise from choosing a different gauge a + 1 
where instead A and B would depend on the gauge parameter. 

Solving the effective Einstein equations for a static source y M = (1,0,0,0) in the Newtonian limit one finds 
that the 00-component of the perturbation h is given by 


- 2 GM 

h oo -- 


43Gjv 

307rr 2 


5G f A(A)-2ff(&) 
127rr 2 r 2 


(4.9) 


which depends on the parameterisation through A and B. The second and third terms in (14.91) are the 
quantum corrections for 0i = 0 due to gravitons and the heavy particle respectively. The crucial step is to 
now find the quantum corrections to the geodesic equation for the test particle of mass m. These arise since 
the particle couples to the quantum metric itself and thus the effective action will also include an extra 
contribution T))) proportional to m. This term is given by 


pnl _ 


/ 


d 4 x^/^g 
A(fa) 


[ 32t r 2 


^m wa log(-n)M^ - ^m^log(-n) [ p ^ + 


^log(-D)Tr + 


B(Pi) 


Rg.M-a)T^ 


(4.10) 


with m^vp C , equal to M M „ pcr with the replacement M -* m and y^ -»■ * v 

In the non-relativistic limit the quantum corrected geodesic equation is given by 

d 2 z 1 1 <5T^ 

d^2 V/l00= mlF 


(4.11) 


where both T^J and /loo depend on the parameterisation. However taking into account both of these 

>2 - 

contributions one finds that -tA is independent of the parameterisation. Explicitly one finds that the 
potential (14.21) is given by 


V(r) 


GnM 

r 


1 + 


43Ga7 

307rr 2 


5 G n 
127t r 2 


7 G n - 

127Tr 2 . 


(4.12) 


where the effect of coupling the test particle is to cancel the parameterisation dependence and contribute 
the last term. Thus we observe that, unlike the beta functions, the Newtonian potential is independent of 
the parameterisation of the quantum fields. 


V. SEMI-CLASSICAL FLOW EQUATION 


We now return to the main topic of this paper namely obtaining a gauge independent beta function for 
Newton’s constant. To this end we will use the results of sections HD and Hill utilising the parameterisation 
(13.51) . The starting point is the functional integral (12.271) on an arbitrary Einstein background with the 
constant mode (f >o removed (in d > 2 this mode will not renormalise the Newton’s constant but should 
modify curvature squared correctio ns in d = 4). Here we shall use the modern formalism of the non- 


perturbative renormalisation group |32-|34|. even though it is not technically needed since we are working 
at one-loop. This approach can then be compared to the old fashioned renormalisation scheme of simply 
evaluating (13.81) and regulating the one loop integrals. To this end we add a regulator term, 


ASfc[<p] = ^ip-TZk ■<£>, 


( 5 . 1 ) 
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to the bare action. Here lZ k is an infra red regulator which suppresses low momentum modes p 2 < k 2 while 
vanishing for high momentum modes p 2 » k 2 . The scale k plays the role of p as the scale which defines 
the beta functions. In the case where we use (13.811 the fields p = {/i 1 ,^, c, c) are all fields appearing in 
the bare action after exponentiating the determinants in functional measure. Alternatively we can use the 
effective actions (13.161) or (13.201) for the unconstrained fields in which case we have p = {h^, a, C M } or 
p = {h^„, C^, Cfj,}. In any case p can be thought of as a vector in field space with components p a and 1Z k 
as a matrix with components 7 Z a b- 

For the covariant momentum p 2 we take the eigenvalues of the differential operators of the operators A„ 
for n = {1,2, s} (given in equations (12.131) . (12.231) and (|3.13l) ~). The regulator then takes a form such that 

S^ 2 \A n ) + TZ k (A n ) = S { 2 \A n ^P n ), P„ = A n + k 2 C(A n /k 2 ), (5.2) 


where C(z) is a dimensionless scalar cutoff function which vanishes for z » 1 but stays finite and non-zero 
for z -* 0. This constitutes a type-II cutoff, in the nomenclature of [35j], where also the potential terms of A n 
are included in the regulator. Another choice, termed type I, is to replace A n by -V 2 in (15.21) . however we 
expect a non-trivial cancellation to occur between the traces of the transverse traceless fluctuations and the 
Jacobian J\ such that only ( d - 3)d/2 propagating degrees of freedom remain. Since, without the regulator 
1Z k , these cancellations must be between traces depending only on the operators A ra , that appear in (13.81) . 
it seems natural that the regulator should also depend solely on these operators as well. In addition it has 
been shown [36j that only type II regulators for fermions correctly regulate modes for the Dirac operator 
in curved spacetime. Finally we point out that it is the properties of the Lichnerowicz Laplacian A 2 that 
determine the stability of classical solutions to the Einstein equations and hence one should expect that this 
continues to be the case in the semiclassical theory. 

After inserting the regulator into the path integral the flowing effective action T k is defined by, 


e -r fc [A] 



e -S[tp]-AS k [cp-(p]+(ip-'fi)-^- 


(5.3) 


where p = (p). One observes that if TZ k diverges for k 2 » p 2 the regulator term becomes a Gaussian peaked 
around p = p, whereas for k 2 -»■ 0 the regulator must vanish and (15.31) assumes the form of the effective 
action. It follows that T^ interpolates between the bare action S for large k and the effective F in the limit 
k -*■ 0. In the semi-classical approximation the flowing effective action takes the form, 


lA - STa] = ^STrlog (S' (2) + IZk) ■ (5.4) 

Due to the relationships between the traces (13.171) . (13.181) and (13.191) and the choice of the regulator (15.21) 
the equivalence between the effective actions T also holds for the k dependent actions 

Tfc - = -Tr 2 r 2 [log (A 2 + 7?.fc(A 2 ))] - — Tiq^ [log (Ai + 7?.fc(Ai))] 

= iTr 2 [log(A 2+ 7e fc (A 2 ))]-Tr 1 [log(A 1 +7e fc (A 1 ))] , (5.5) 


where we have absorbed irrelevant constant factors of 167 tG{, into the fields. Setting lZ k = 0 we recover 
Taking a derivative with respect to the Wilsionian RG time t = log/c/fco to obtain the following gauge 
independent flow equation, 


o r _1 T , d t TZk{ A 2 ) 

‘ k 2 Tl2T * A 2 +lZ k (A 2 ) 


1 JMA^ 

2 u A 1 +TZ k (A 1 ) 


d t TZ k ( A 2 ) 


= - IV ■ 

2 A 2 + 7Z k (A 2 ) 


- Tr , 


dtTZkiAj) 
Ai + 7Z k (Ai) 


(5.6) 


Note that the flow equation (15.61) . unlike (13.81) . is both UV and IR finite and therefore does not necessitate 
the introduction of an explicit UV cutoff. Instead the bare action enters the flow equation as a boundary 
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condition. Including the constant mode (f >o we give an extra term in he first line of (15.61) arising from this 
single mode. In the second line of (15.61) we give the form of the flow equation for unconstrained fields. 

Here we consider the following ansatz the flowing action 


T k = 


1 


167rGfc 


J d d xy/g(2Xk - R) = S[gp V ] + J d d x%/g (5Xk - dnuR) , 


(5.7) 


which takes the Einstein Hilbert form and allows us to extract the beta functions for the running couplings 
Gk and A k- Equivalently we define the vacuum energy density A k and the inverse Newtons couplings Kk by 


A k = 


A k 


X b 


8nGk 8nGb 


■ dX k 


Kk = 


1 


1 


167rGfc 167rGb 


- 5k 


k ■ 


(5.8) 


The UV scale A can then be defined as the scale for which SXk=A = 0 = Snk=A- In turn these boundary 
definitions can be viewed as replacing the operational meaning of the ‘bare quantities’ Gf, and At, entering 
the functional integral. Thus, the flow equation defines a renormalisation scheme whereby the UV boundary 
condition replaces the bare action and the IR boundary condition sets the renormalisation condition in terms 
of renormalised quantities obtained in the limit k = 0. The advantage of this scheme is that it dispenses 
with with the formally divergent path integral and generalises beyond the perturbative regime 32 34]. 

We now put the flow equation together expressing all quantities in units of k in order to find the au¬ 
tonomous system of beta functions. In particular we drop the index k from the couplings when referring 
to the dimensionless quantities such that G = k d ~ 2 Gk and A = k~ 2 Xk- The dimensionless flow equation then 
takes the form 




d t G- (d- 2)G 
16t rG 2 


R + d t X + dX = Tr 2 [W(A 2 )] - 2Tn[W(Ai)]. 


where we define the function, 


W(z) = 


C(z)-zC'(z ) 


z + C{z) 

This gauge independent flow equation (15.91) is the main result of this section. 


(5.9) 


(5.10) 


VI. HEAT KERNELS AND UNIVERSALITY 

To obtain the beta functions one computes the traces on the RHS of the flow equation (15.91) and compares 
the terms on each side of the equation to linear order in the curvature. In the same way one way evaluate 
(13.81) directly to determine the counter terms needed for the one loop effective action on an Einstein space. To 
evaluate the traces in both cases we therefore exploit the early time heat kernel expansion for the operators 
Ai and A 2 acting on the vectors and symmetric tensors respectively. For a general function /(A n ) one has 
the following expressions, 

Tr„[/(A„)] = —— -Y,Q±-iU] / d d Xy/ga 2i ,n, (6.1) 

(4tt)2 i=Q 2 J 

where the functionals Q m for m > 0 are given by 

Qm[f] = -^ Hdzz^fiz), ( 6 . 2 ) 

1 (m) Jo 

and the coefficients A 2 i in are the heat kernel coefficients for the operators A n proportional to curvature 
monomials with mass dimension 2 i. In order to find the beta functions for G and A we need only the heat 
kernels for i- 0,1. For i = 0 the heat kernel coefficients are equal to the dimension of the field space 

d(d+ 1) 


ao,i - d , ao, 2 


2 


(6.3) 
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whereas for i = 1 the heat kernel coefficients of a generic operator of the form 

A„ = -(V 2 + E n ), with (E 2 )^^ = 2R^ q ^, (E 1 r a = R% (6.4) 


are given by 

a 2 ,n — ~zd n R + ti n E n (6.5) 

6 

where the trace tr„ denotes the trace over the field space leading to tiiEi = R and tr 2 £ , 2 = -R We then 
obtain 

Tr 2 [W(A 2 )] - 2Tr[W(A!)] = ((a 0 ,i - 2a 0 , 2 ) Q- [W] + (^N g + tr 2 E 2 - 2^^ Q^W]} (6.6) 

= 7 -^ [ d d x^/g(N g Q AW] + \(N g -18) Q, 1 [W]r) + ... (6.7) 

(47t)2 J \ 2 D 2 / 

where iVg = ^d(d- 3) is the number of propagating degrees of freedom in d-dimensional quantum gravity. 
One sees here the universal origin of the factor N g -18 coming just from the heat kernel coefficient. We note 
that one could of have simply renormalised via the introduction of counter terms without the introduction 
of the IR regulator. That is we could evaluate (13.81) directly to obtain 

F = S+ 1 f d d x^/g[N g Q A\og] + UN g -18) Qd Jlogji?) +... (6.8) 

2(4-7t)2 J \ 2 o 2 / 

Then one simply observes that the counter term for the v /yi?( 7 ) term in (12.21) would be proportional to 
N g -18 and the volume term is proportional to N g . Thus we see, from this point of view, why it is important 
to make the regulator a function of the differential operators appearing in the hessians themselves. Any 
other choice, such as 7 Zk = 7?.fc(-V 2 ), would lead to non-universal heat kernel coefficients that do not appear 
in the one-loop effective action. One should also stress that the origin of -11 +1 Nf can also be traced to the 
heat kernel coefficient for the corresponding tensor structure in Yang-Mills [37j. In addition one my also go 
to the next order in the curvature expansion where in d = 4 the result is universal since Qo[W] = W(0) = 1. 
Here we find explicitly the curvature squared term 

Tr 2 [W(A 2 )]-2Tr 1 [W(A 1 )] = ... + ^ip f d d x^g(^R^R^-^R^ + ... (6.9) 

reproducing the one-loop counter term found in [2(J but with the cosmological constant replaced by A& -*■ 
R/4. Here we note that the inclusion of the constant mode contribution (which does not contribute the the 
renormalisation of Newton’s constant) is need to correctly evaluate these terms. 

We note the general pattern that universal factors appearing in beta functions and counter terms through¬ 
out quantum field theory can be found by use of heat kernel techniques. These factors can be found for 
dimensionful couplings, such as Newton’s constant, as well as dimensionless ones. The difference is that the 
renormalisation of dimensionless couplings do not come additionally with non-universal factors. However 
the sign of non-universal factors appearing in the renormalisation of dimensionful should be universal and 
gauge independent if they are to be physically meaningful. The fact that the universal factor found here 
only depends on the dimension through the number of polarisations would seem to support the view that 
this is indeed the case. 

It should be noted, however, that the unique factors N g - 18 and N g could be modified if one was 
to regulate the auxiliary fields fields differently than the transverse traceless fields. Thats is if we choose 
different regulators Ci(z) and C^C-s) for the separate traces in (15.91) we would end up with different functions 
W(z). Making such a choice would mean that the factors N g and N g -18 of (16.61) no longer be forthcoming. 
Equally if we choose to regulate with a type I cut off such that 7 Zk = 7 Zk(-V 2 ) rather than the choice (15.21) 
these factors would again not appear in (16.61) . One must conclude then that if the factor N g - 18 is truly 
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universal and physical only a sub set of all possible regulator choices leads to physical results. What we 
suspect is that regulator schemes must break diffeomorphism invariance sufficiently mildly such that the 
cancellations of unphysical polarisations continues to be present. 

To back up our point of view imagine we pick C\(z) + C 2 (z), then the term in (16.61) proportional to the 
space time volume would not necessarily be proportional to d- 3. This being the case the vacuum energy 
would appear to be renormalised in d = 3 even though the evaluation of the unregulated effective action 
(16.81) would lead to no divergencies proportional to the volume. A similar argument can be made for the 
factor N g - 18 in the critical dimension d = d c for which N g (d c ) = 18. Thus only by picking the regulator 
scheme as we have do we avoid such inconsistencies. 


VII. BETA FUNCTION FOR NEWTON’S CONSTANT 


It is now straightforward to find the beta function (3g = dtG for Newton’s constant by inserting (16.61) into 
(15.91) . For general cutoff function C(z ) it reads 

fe - (d - 2)G+ l (4^uV) ^ l[C1G2 (7 ' 1> 

where the regulator dependent functional X n [C] is given by the integral 


T 

J -'Tl 


[Ck] = f 

Jo 


dzz n ~ 1 W(z) = f°°dzz d/2 - 2 
Jo 


C(z)-zC'(z) 
z + C(z ) 


(7.2) 


The first term in (EH) is the classical scaling arising from the dimensionful nature of Gk ■ The second term 
is a quantum correction proportional to the momentum integral 2 d/ 2 -i[C] indicating that this equation is 
one loop. 


A. d = 2 + e and d = 2 dimensions, a tale of three beta functions 

Before continuing to the general d case we take a digression into the behaviour of the beta function near 
two dimensions and try to understand the subtleties of the limit d -*■ 2. First we note that in the limit d -*■ 2 
we have 

2 d /2- i [ C ']| d _ 2 = ^2 + finiteterms ( 7 - 3 ) 

where the singular term is independent of the regulator function C. However the coefficient of G 2 has a 
finite universal limit for d -»■ 2 owing to the presence of the gamma function. This allows one to extract a 
universal beta function in d = 2 + e dimensions and in the case d = 2. 


1. d = 2 + e dimensions 


For d = 2 + e we find the universal beta function 




(7.4) 


We note that this is in agreement with previous studies 38-4c| using the linear parameterisation (12.51) but 
differs from the result obtained with the exponential parameterisation (12.61) jfl 41] with p = d/2 - 1 and 
\b = 0 which gives 


P G = eG-|G 2 . 


( 7 . 5 ) 
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The reason for this is that by working off-shell and setting Ab = 0 the exponential parameterisation gives a 
hessian for the scalar mode (j> becomes oc (d- 2) A and thus the cancellation between these modes and 
Jo does not occur. Instead the hessian for the conformal modes mimics the induced action 


& 


f d 2 x^/g(jR + 


16 


(7.6) 


of the conformal anomaly in d = 2 despite being of order e [41]. In the equations leading to (17.41) such 


fluctuations are not present since all conformal fluctuations are cancelled by the functional measure. The 
result m can then be understood as arising from the particular way the limit d -*■ 2 is taken where the 
fact that the hessian for <j> vanishes is not taken into account. The limit then paradoxically reproduces the 
expected result in d = 2 in the presence of the conformal anomaly even though the induced action (17.61) for 
has not been accounted for. 

In order to help clarify the situation we briefly recall the derivations of m given in |id| and generalise 
it to an arbitrary parameterisation. The derivation differs in two respects from the methods used in this 
paper. Firstly dimensional regularisation is used such that only the logarithmic divergences are retained 
and secondly the freedom to redefine the background metric is exploited in order to define a physical beta 
function. Thus first one obtains the one loop equation 


dd'[g» v ] = 




(7.7) 


where b = b(a,p,/3i) is a number which depends on the gauge (12.91) and parameterisation (14.31) and /i is the 
renormalisation scale. For the special case (13.51) one has 6 = 0. To remove the unphysical dependencies in 
the case b * 0 the background metric is redefined g -»■ g' = Z g g^ in such a way that the cosmological 
constant with respect to g' is not renornralised. Explicitly one has 


Z~ 1 dtZ g = &KGffeb 


(7.8) 


leading to 


dtT[g^] = p e J d 4 xy/g J 


16-7T 3 


(7.9) 


reproducing JO) independently of both the gauge and the parameterisation. At a deeper level one finds 
that the renormalisation of the dimensionless product G A ^+2 is independent of Z g . Indeed one may define 
the beta function via the relation between the bare and renormalised dimensionless products in units p - 1 


GK- 2 = ffV- 2 (i + fA_o) 


(7.10) 


which is independent of b and hence both the gauge and parameterisation. The beta function is then found 


by asking how G should be changed if A A+ 5A while keeping G^X^ 2 fixed according to the definition [401 


SG = -~0g- (7.11) 

a X 

This beta function is again given by m - Thus we observe that when considering the physical renormali¬ 
sation there is no dependence on the parameterisation. The advantage of using a parameterisation or gauge 
for which 6 = 0 is that the naive beta function automatically agrees with the physical one. One should note 
that different physical beta functions may be defined once other couplings or masses are introduced [4o| . 


2. d = 2 dimensions 

One should note that the beta function (17.41) is valid for d + 2 and that for the strictly d - 2 the hessian 
for the modes <j> vanishes reflecting conformal invariance. To work explicitly in two dimensions we note 
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that f drXy/jR is a topological invariant however the inclusion of both gauge fixing and ghost terms means 
that the path integral is non-trivial. Furthermore since there are no transverse-traceless fluctuations in two 
dimensions the use of the transverse-traceless decomposition is problematic. As such we will use de Donder 
gauge (i.e. for a = 1 and p = 0) such that the hessians and ghosts determinants are given (13.131) . (13.221) and 
(13451) . 

We observe then that there are both traceless metric and ghost contributions while the hessian of the 
conformal fluctuation vanishes. Since the second variation of the Einstein-Hilbert action itself vanishes 
in d = 2 it is the case that (13.221) is the second variation of the gauge fixing action alone even though in 
higher dimensions it results from a cancellation of gauge fixing terms with terms in the second variation of 
Seh = -f d 2 x^R. To make sense of this we can observe that an explicit computation of = 0 for d = 2 
gives a relation which implies that the second variation of the gauge fixing action is 

16tt G b h-S™ hh -h=- f (Pxy/gh^VpVxh^lf d 2 x^gh^ , (7.12) 


leading to (13.221) . As a crosscheck the second equality above can be confirmed by writing h flL , = V 




- g /ji/V a v a which is a complete decomposition of a traceless symmetric two tensor in two dimensions. 
If we additionally include D scalar fields these can be interpreted as the D spacetime dimensions in string 
theory. That is we include an action of the form 


String = f d d X ^Gud^d^X" 


d-> 2 


(7.13) 


with I = where we assume G/j to be a flat metric over the D dimensional spacetime in which the 

string lives. Going through the motions one then obtains a beta function for d -*■ 2 given by 


0G ~ -(26 - D) • ^G 2 . 


(7.14) 


The difference between (17.41) and (17.141) arises from the lack of conformal fluctuations contributing to the 
later which therefore do not cancel the Jacobian Jo which in turn leads to an additional term in the flow 
equation along with a contribution from the D scalar fields. The beta function (17.141) is the expected result 
reflecting the fact that when 26 scalar fields are added to the action the conformal anomaly vanishes. Thus 
we reproduce the result that the bosonic string lives in D - 26 dimensions, in which case we do not need 
to bother about the conformal anomaly [42j]. If D < 26 the beta function (17.141) does not vanish and the 
conformally anomaly should be taken into account by including the induced action (17.61) . This has been 
investigated recently in [43j where it was found that in d = 2 the ghosts alone gives the “-26” contribution 
found whereas here the “-26” arises from a combination of ghosts and gauge fixing. This is nonetheless 
consistent with our result (17.141) since the conformal fluctuations which are present there essentially play 
the role of one of the D scalar fields. 

We conclude that (17.41) is the physically meaningful gauge independent beta function for e > 0 in agreement 
with the physical beta function defined by (17.111) . On the otherhand in strictly d = 2 one needs to take 
into account the conformal invariance (or anomalous breaking thereof) to correctly evaluate the functional 
integral and find the corresponding beta function (17.141) . The difference between (17.41) and (17.51) arises since 
dm is the result of quantising general relativity in d > 2 where conformal fluctuations do not propagate 
and are exactly canceled by the Jacobian Jo. The beta function (| 7.5 1) arises when working in the limit d -> 2 
where, due to the particular choice of gauge/parameterisation, the conformal fluctuations mimic the effect 
of the induced action (USD such that the cancelation with Jo is no longer exact. 


B. d-dimensional beta function for Newton’s constant 

In d > 2 dimensions the integral Td/ 2-1 [G] depends on the regulator function. This must be the case since 
Newton’s coupling G& is dimensionful and thus G will always depend on the regulator at least up to rescalings 
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of k. However it is interesting to note that for an optimised cutoff [44[ of the form Copt = (zq ■ 
where ©(a;) is a Heaviside theta function and zo is a positive constant, one obtains 

2 


■^d/2—1 [Copt] _ Z 0 


d -2 


z)Q(z 0 -z), 


(7.15) 


valid in all dimensions d > 2. Hence a property of the optimised cutoff is to set all higher order corrections 

d -2 d 

in (17.31) to zero up to an arbitrary rescaling. For z 0 2 = (47r)2^ 1 r(^) one then obtains the beta function 

/3 g = (d-2)G-^(18-N g ) G 2 , (7.16) 

where only the universal factor of b remains. This form of this beta function is clearly universal up to the 
normalisation of the RG scale k or equivalently G. Thus although G itself is not a physical observable it 
is clear that the sign of the quantum correction to /3g is universal since one has Z„[G] > 0 for all regulator 
functions C. 

For N g < 18 the beta function implies Gk decreases as the scale k is increased. The running of G stops at 
a UV fixed point given by, 


d-2 3(4 7 r)^ 1 F(^) l 
2 18 - N g I d / 2 -i [C]‘ 


(7.17) 


This fixed point exists for positive G* in all dimensions for which N g < 18 and describes an asymptotically 
safe quantum field theory at high energies. However since this is a semiclassical result this interpretation 
is rather premature and one would like to go beyond this approximation to confirm this conclusion. Indeed 
the critical exponent 

l/v=-^(G*) = d-2, (7.18) 


which describes how G approaches G* is simply the canonical one, independent of the regulator. To compute 
quantum corrections to v one must go beyond the semi-classical approximation by computing /3g to higher 
orders in G. This will be investigated in a companion paper 45[ where we exploit a non-perturbative 
approximation to find quantum corrections to v. 


C. Paramagnetic dominance 

From the form of the beta function (17.161) we observe that the increase of spin two degrees of freedom is 
responsible for the loss of asymptotic safety in higher dimensions. The term -18 which allows for asymptotic 
safety in d < 8 dimensions is universal and independent of the dimension. The origin of this term is the 
interactions of the metric fluctuations with the curvature of spacetime, and represents the ‘non-Abelian’ 
nature of gravity, in agreement with the ideas put forward in [dcj. In particular it’s origin is the ’non- 
minimal’ coupling terms E n of the differential operators (16.41) 

Following |46[ we then denote terms proportional to the non-minimal couplings E n as ‘paramagnetic’ and 
those arising from the Laplacian as ‘diamagnetic’ to obtain the expression 

A2, n -2A2,i = (d(d + l)/2-2d-6(tr 2 E2-2tT 1 E 1 )) ^ (7.19) 

\_ v _✓ _ v _/ 6 

diamagnetic paramagnetic 

= ( ^ ) f • (7-20) 

i. paramagnetic 

diamagnetic ^ ° 

This confirms the results of [46[ in a gauge independent setting: the physical mechanism behind asymptotic 
safety is related to the metric fluctuation’s paramagnetic interaction with the curvature of spacetime. This 
effect is countered in d > 3 dimensions by the diamagnetic interactions encoded in the Laplacians which 
dominate starting in d > 8 dimensions. 
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VIII. RENORMALISATION OF THE VACUUM ENERGY 


As we have seen the absence of the cosmological constant in the beta function for Newton’s constant 
results from demanding gauge independence and therefore its presence in gauge dependent beta functions 
is most likely unphysical. This can be seen as resulting from the fact that the cosmological constant is not 
a mass for the graviton. Rather the vacuum energy only couples to the conformal fluctuations a which are 
not themselves dynamical. Nonetheless the vacuum energy A = has also a k dependence given by, 


(){A — dX + 


N g 

(47r)*r(d/2) 


1AC] 

2 


( 8 . 1 ) 


The quantum correction being proportional to N g . This beta function has a fixed point, 


A* 


1 

d 


N 

- T 1 -[G] . 

(47r)2T(d/2) 2 


( 8 . 2 ) 


Note that in the absence of any propagating degrees of freedom, i.e. when taking d = 3, the quantum 
corrections to this beta function vanish. This result is just related to the divergencies of the vacuum 
encountered in any quantum field theory in flat spacetime. In this case we are just looking at the graviton 
contribution to the vacuum energy. To see this lets write (110.41) explicitly in terms of the IR cut off scale k 
and the dimensionful vacuum energy A& = k d A, then one has 


d t X k = dk d X* , 


(8.3) 


introducing the bare UV scale A defined by Xk=A = Ab we can solve this equation to find, 

Xk = {k d - A d ) A* + A b . (8.4) 

In the limit k -> 0 we remove the IR regulator to obtain the observed vacuum energy Ao = Afe=o- In turn we 
may express the bare vacuum energy in terms of Ao as 

Ab = Ao + A d A*. (8-5) 

Thus the ‘beta function’ (110.41) is nothing but the statement that we must include a counter term proportional 
to the number of local degrees of freedom in order to get a finite renormalised vacuum energy Aq. 


IX. CONTINUUM LIMIT 


As discussed previously the existence of a continuum limit, A -* oo, relies on the sign of b , the G 2 coefficient 
of (3g- To understand this we note that the sign of the bare Newtons constant Gb must be positive if the 
functional integral is to make sense. If one is forced to take Gb < 0 to renormalise the theory one must 
concede that the theory cannot be fundamental and that the continuum limit does not exist. From the (ED 
one may infer the energy dependence of the Planck mass Kk = ; 

Kk = (k d - 2 -A d ~ 2 )K* + K b , (9.1) 


where k* = 16 J G is the regulator dependent value of the fixed point for the dimensionless quantity n - 
k 2 ~ d Kk- Taking k -* 0 one has the renormalised Planck mass Mp^ 2 = — = Kk=0 where Gjv is the measured 
value of Newtons constant. Then we have 


1 1 


167rGb 167rGjv 


+ A d 2 k* . 


(9.2) 













23 


for the functional integral to make sense one must have Gb > 0 otherwise the kinetic terms will have the 
wrong sign and the theory will breakdown. In particular for the semiclassical theory to be valid we need 


1 

167tGjv 


> -A 


d- 2, 


(9.3) 


Now we know that tc* oc (18 - N g ) and that the proportionality constant is positive from which we obtain 
the inequality for N g > 18 


/ A \ d 2 

iwJ <16 ’ r|G - 1 ' (9 - 4) 

Thus, while for 18 < N g there is no maximum value for which the UV scale can take, for N g > 18 (USD 
provides a bound. In the latter case we can say that the semiclassical theory predicts its own downfall at 
a finite energy scale. In the former case no bound exists and the semiclassical approximation suggests that 
the theory is asymptotically safe. 


X. INCLUSION OF MATTER 


It is rather straight forward to include the effects of matter. This is done by adding a matter action to 
the bare Lagrangian and adding a corresponding regulator term (see for example [13] ) or by computing the 
regulated one loop effective action 48|. Doing the former adds new traces to the right hand side of (15.61) 
proportional to the number of matter fields. In turn these give contributions to the beta functions. Here 
we give the result in d = 4 spacetime dimensions, although the result is easily generalised to arbitrary d , 4 


Pg = 2G - -(16 - 2N d + AN M - N S )G 2 , 


( 10 . 1 ) 


up to the normalisation factor, where Np is the number of Dirac fermions, Nm is the number of gauge fields 
and N s is the number of scalars. Plugging in the values for the standard model (Np = 45/2, Nm = 12, N s = 4) 
we have 


Pg = (d-2)G - ^ • 15 G 2 , 


( 10 . 2 ) 


which continues to be anti-screening as in the pure gravity case. We therefore can conclude that the semi¬ 
classical theory of the standard model coupled to gravity is not predicting its own down fall at the quantum 
level and the existence of a continuum limit along the lines of asymptotic safety remains a possibility. In 
particular the beta function predicts asymptotic safety provided 


16 — 2 Np + 4Nm ~ N s > 0 


(10.3) 


otherwise G* is negative and the bound (19.41) applies. We observe that by adding large numbers of fermions 
or scalars the theory can break down a fraction of the Planck scale. This then puts constraints on effective 
theories that contain a large number of scalars or fermions. 

Similarly the renormalisation of the vacuum energy at one-loop is given by 


dt A — —dX + 


-‘’Bose AVerini 

(47r)Tr(d/2) 


Id [G], 


(10.4) 


where the number of bosons Nb OS e (including gravitons) and the number of fermions fVFermi enter with 
opposite signs. This is the expected result and entails that the vacuum energy is not renormalised at 
one-loop for super-symmetric theories. 


4 Here we assume the fields are massless and minimally coupled to gravity. The beta function will be modified if these 
assumptions are dropped. 
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XI. CONCLUSIONS 


In this paper we have derived a gauge independent effective action for quantum gravity at one-loop which 
only depends on physical fluctuations of the metric. To achieve gauge independence we have used a specific 
parameterisation of the metric for which the volume element is linear in the conformal fluctuation a. We have 
then used renormalisation group techniques to compute the one-loop beta function for Newton’s constant. 
When introducing a regulator scheme we have been careful to do so in such a manner that cancellations, 
which should occur when diffeomorphism invariance is present, continue to occur in the presence of the 
regulator. In particular this relies on choosing a type II regulator of the form (15.21) while choosing the 
same regulator function C(z) for the different fields. Different choices can lead to unphysical results. In 
particular the absence of the local fluctuations in d = 3 may not be manifest if the regulator functions are 
chosen differently. 

Here we offer an interpretation of the gauge independence of our results. A key point is that our calcu¬ 
lations are performed on Einstein spaces for which the trace-free Einstein equations are solved. Thus we 
are on-shell with respect to all but one of the equations of motion, this one being essentially the equation 
of motion for the volume element. Now, one way to ensure gauge independence in any gauge theory is to 
work in terms of physical quantities rather than gauge variant fields such as the metric 49j. Since we single 
out the volume element by our parameterisation we are therefore singling out the corresponding physical 
quantity which in this case is the gauge invariant field cf> given in (13.71) . Thus we solve the problem of finding 
a physical parameterisation on an Einstein space (at least at one-loop order). This is seen explicitly since 
without the gauge fixing terms the quadratic action involving S^ is independent of both £ and ip and hence 
diffeomorphism invariant. 

A further insight is gained by noting that the effective field equations for the field a are solved for the 
expectation value (a) which then gives direct access to the mean space-time volume 


(V('y)) = V(g) +f d^xy/g^-. ( 11 . 1 ) 

This can be considered as advantage since to compute (V('y)) for the standard effective action involves an 
infinite number of terms consisting of n-point functions to all orders by expanding the volume element ^7 
in h^. We therefore believe that the general scheme put forward here is physically well motivated and 
should be used in approximations beyond the simple one employed here. The fact that the universal factor 
m only depends on the dimensionality of spacetime through the number polarisations N g would seem to 
be a vindication of this claim and that it represents the analogy to the factor 11 - |Nf in QCD. 
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